The criterion for dividing the surrounding rock EDZs of underground caverns based on the energy dissipation degree

An energy calculation parameter named the energy dissipation degree (RUd) is introduced based on the analysis of the energy dissipation mechanism and energy evolution characteristics during conventional triaxial tests of the granite of Shuangjiangkou. The deviatoric stress‒strain curve of rock can be divided into five stages using four stress thresholds (crack closure stress σcc, crack initiation stress σci, damage stress σcd and peak stress σp), which also correspond to the four RUd thresholds (RUdc, RUdi, RUdd and RUdp) on the energy dissipation degree–strain curve. A given stress threshold increases with increasing confining pressure; however, a given RUd threshold is basically stable under different confining pressures. Then, a new criterion for dividing the excavation damaged zones (EDZs) in the rock surrounding underground caverns based on the monotonically increasing characteristics of the energy dissipation degree‒axial strain relationship curve is proposed, and it allows for the classification of the surrounding rock into five types of zones through quantitative analysis of the RUd thresholds. Based on the criterion for dividing the EDZs of the surrounding rock mass of the underground cavern, the EDZs of the surrounding rocks of the underground cavern group of the Shuangjiangkou Hydropower Station are analyzed. The distribution characteristics of the EDZs of the rock surrounding underground caverns obtained by numerical simulation calculations based on RUd are basically the same as those obtained by in situ elastic wave tests. However, the RUd-based method for classifying the EDZs of the surrounding rock has the obvious advantage of being able to probe the boundaries of the undamaged zone (UDZ) of the surrounding rock more explicitly, while the method based on wave velocity testing is not sufficiently explicit. The damage zoning of the surrounding rock based on RUd can provide support design advice for the excavation of the surrounding rock, such as the support method, the length of the free section and anchor section of the prestressing anchor, etc.


Introduction
During the construction of underground projects, rock excavation and the protection of surrounding rock mass stability are conflicting but necessary challenges. When using static and dynamic rock breakage methods (soundless cracking agents, tunnel boring machine, gas blasting, rock blasting, etc.) to excavate rock, the excavation and unloading process will disrupt the ground stress balance and cause redistribution of stress in the surrounding rock [1][2][3][4]. The energy stored in the surrounding rock is gradually dissipated, while geological, thermal, mechanical, hydraulic, chemical and acoustic characterization of the surrounding rock are affected [5][6][7]. The damage and destruction of the surrounding rock inevitably occurs, resulting in excavation-disturbed zones or excavation-damaged zones (EDZs) [8]. The temporal and spatial evolution of EDZs are characterized by the gradual development of the surrounding rock from a relatively intact state (in situ stress) to a fully adjusted state (new stress equilibrium state after excavation) under the influence of excavation [9]. The formation and expansion of EDZs are the root cause of deformation and the damage of the underground cavity envelope, causing challenges related to stability and seepage, thus impairing the performance and function of the surrounding rock mass [10]. Therefore, it is critical to select reasonable and practical methods to accurately determine the formation mechanism, change trend, and the range of the EDZs of the surrounding rock to provide support for the design of underground cavern support.
To better conduct research on EDZs, many countries or regions have set up specialized underground research laboratories (URLs), such as the URL for the Lac du Bonnet granite in Canada [8,9], the Mont Terri Rock Laboratory for the Opalinus Clay in Switzerland [10,11], the Ä spö Hard Rock Laboratory for granite in Sweden [12], the URL for the Méraillet crystalline rock in the French Alps [13], the Tournemire URL in Francev [14], the Meuse/Haute-Marne URL in France [15][16][17], the High-Activity Disposal Experimental Site in Belgium [18], the Tono mine in Japan [19], the Horonobe URL in Japan [20,21], the KAERI Underground Research Tunnel in Korea [22] and the Jinping URL in China [23,24].
To quantitatively determine the range and variation pattern of EDZs, scholars have conducted studies based on various instruments to record relevant data, such as displacement and strain measurement [19,25], wave velocity testing [24,26], fluid permeation testing [21,27,28], gas permeation testing [5,29], resistivity testing [12,30,31], acoustic emission and microseismic event monitoring [32][33][34], borehole core analysis [13,35], chemical test methods [36,37], seismic wave testing [38][39][40][41], borehole camera technology [42,43], drill hole panoramic digital imaging technology [24], filling resin observation [17,44], and physical model tests [33]. It has been shown that the characteristics of EDZs vary with the rock properties, in situ stress field, geological environment, excavation method, and cavity section geometry, and the results obtained from different research methods vary. The relevant data recorded by the instruments reflect the unique thermal, mechanical, hydrodynamic, chemical, and acoustic characteristics of EDZs, which can be used to gain insight into the mechanisms of EDZ formation, and to determine the extent and variation patterns of EDZs. Nevertheless, owing to topographical and financial constraints, among other factors, the amount of data recorded by instrumentation is limited, challenging the assessment of the extent and degree of EDZs. Numerical simulation methods have also been widely used for the study of EDZs [45][46][47][48]. In addition, some scholars use reliability analysis tools to estimate the probability of damage zone generation for the optimal design of support measures [49]. These methods for generating quantitative results are often used in combination to improve the accuracy of the methodology of the study [12,14,22,25,48].
Based on the analysis of the energy variation in the Shuangjiangkou granite via conventional triaxial tests, the energy calculation parameter R Ud is introduced in the paper. Then, the criterion for dividing the EDZs in the rock surrounding underground caverns based on R Ud , which considers the effects of both confining pressure and strain, is proposed. The surrounding rock EDZs of the underground caverns of the Shuangjiangkou Hydropower Station are classified by applying the criterion through numerical simulation. In the 0+070.00 m section, eight in situ wave velocity test boreholes have been arranged in the main powerhouse, the main transformer chamber and the tailwater surge chamber, and the EDZs of the surrounding rock mass are classified according to the variation in wave velocity with the depth of the borehole. The EDZs based on elasticity wave testing data from in situ tests are compared with the results obtained from numerical simulations. Analyzing the surrounding rock's overall stability through zoning can provide a reference for the design of support.

Conventional triaxial tests
Shuangjiangkou Hydropower Station is located in the upper reaches of the Dadu River in China, and the rock in which the underground powerhouse engineering is located consists of granite. The deviatoric stress-strain relationship of granite under different envelope pressures is obtained by conventional triaxial tests, shown in Fig 1 [50].

The energy dissipation degree R Ud
According to the principle of thermodynamics, the total energy received by granite samples during tests can be calculated by Eq (1) [51,52].
where U, U l , U e and U d represent the total absorbed energy density, freezing energy density, elastic strain energy density and the dissipated energy density, respectively. Since the amount of freezing energy is remarkably small compared to the elastic strain energy and the dissipation energy [52], an approximation is taken in the energy calculation, ignoring the effect of freezing energy on the deformation damage of the rock sample, so Eq (1) is simplified to Eq (2), Eq (3) can be used to calculate the total absorbed energy density, where U 0 is the energy density caused by hydrostatic loading, U 1 is the energy density owing to axial compression by σ 1 after hydrostatic pressure, and U 3 is the energy density owing to lateral compression by σ 3 after hydrostatic pressure.
the dissipated energy can be obtained by Eqs (2)-(7) [50,51]. The rock damage evolution process will lead to energy dissipation, and the study of the damage evolution characteristics of rocks from the perspective of energy dissipation can fundamentally reflect the deformation and damage characteristics of rocks. To better study the energy evolution of granite during the damage process of the Shuangjiangkou Hydropower Station, an energy calculation parameter named the energy dissipation degree (R Ud ) is introduced. R Ud is defined as the ratio of the dissipation energy during the evolution of rock damage to the dissipation energy at the peak stress, and R Ud represents the degree of dissipation energy change during the damage of the rock, which can be obtained by Eq (8). A large number of studies [53][54][55] have shown that, the deviatoric stress-strain curve of rock can be divided into five stages using four stress thresholds: crack closure stress (σ cc ), crack initiation stress (σ ci ), damage stress (σ cd ) and peak stress (σ p ). The four stress thresholds correspond to the four R Ud thresholds (R Udc , R Udi , R Udd and R Udp ) on the energy dissipation degreestrain curve, as shown in Fig 3. The determination of the energy dissipation degree thresholds is shown in Table 1 and The stress threshold at different confining pressures becomes larger with increasing confining pressure, while the corresponding threshold of energy dissipation degree R Ud varies in a smaller interval at different confining pressures.

Definition of the criterion for dividing EDZs based on R Ud
In analyzing the EDZs of the surrounding rock encountered in underground cavities, stress, deformation, and plastic zone variation characteristics are generally taken as the foci of research. In contrast, less research has been conducted to analyze EDZs based on energy dissipation characteristics [47,56]. The damage evolution characteristics of rocks are studied from the perspective of energy dissipation, which can reflect the deformation process of rocks.

PLOS ONE
The criterion for dividing the surrounding rock EDZs of underground caverns

PLOS ONE
The criterion for dividing the surrounding rock EDZs of underground caverns In underground cavity excavation, the initial ground stress is unbalanced owing to the excavation activities, which redistribute the stress of the surrounding rock. In general, the stress concentration will occur tangential to the cavern wall, whereas stress relaxation will occur along the cavern wall's radial direction, and even tensile stress will appear. Such unfavorable stress conditions result in varying degrees of rock damage, and the destruction around the cavity chamber and the surrounding rock can be classified as five circles according to the distance from the cavity chamber, from near to far (Fig 5).
The construction-damaged zone (CDZ) is generally located immediately around the underground cavity chamber, where the surrounding rock's damage and destruction is serious, while numerous macroscopic cracks or fracture surfaces are formed, and the surrounding rock loses its bearing capacity. The surrounding rock in this part can be supported by initial protection with shotcrete immediately after excavation, and then supported by system anchors with reinforcement mesh and poured concrete lining. The pressure generated in this area is loose ground pressure and should be used as the basis for the design of the pressure to which the support is subjected. The highly damaged zone (HDZ) is generally located in the second circle around the underground cavity chamber, in which the surrounding rock's damage is serious, and additional macrocracks or rupture surfaces occur, seriously weakening the surrounding rock's bearing capacity. The pressure generated in the area where CDZ and HDZ are located is loose ground pressure and should be used as the basis for the design of the pressure to which the support is subjected. The mildly damaged zone (MDZ) is located in the third circle around the underground cavity, in which the surrounding rock's damage is slight and there are fewer macrocracks or rupture surfaces; thus, the surrounding rock has a certain bearing capacity. The pressure generated in the area where MDZ is located can be used as the design basis for the pressure on the support, and it should be used as the key monitoring part for the stress and displacement of the surrounding rock, and then the working condition of the surrounding rock and support can be monitored in real time. The excavation-influenced zone (EIZ) is located in the fourth circle around the underground cavity, where the surrounding rock is influenced, there are fewer cracks, and the surrounding rock is of better quality, with fewer cracks and stronger bearing capacity. The area where the EIZ is located can be used as a key monitoring site for the surrounding rock stress and displacement, and then the working condition of the surrounding rock and support can be monitored in real time. The undamaged zone (UDZ) includes the very slightly disturbed area and the original rock area. When prestressed anchors are used for support design, the area in which CDZ, HDZ, MDZ and EIZ are located should be designed as the free section of anchors, while the area in which UDZ is located can be designed as the anchored section of anchors.
Based on conventional triaxial tests, the deviatoric stress-axial strain relationship curves and the energy dissipation degree-axial strain relationship curves are used as the basis for dividing the EDZs of the surrounding rock. This is the criterion proposed for identifying the EDZs in the rock surrounding underground caverns based on the energy dissipation degree R Ud , as shown in Fig 5. The criterion uses four energy dissipation degree thresholds (R Udc , R Udi , R Udd and R Udp ) to classify the surrounding rock mass into five categories: CDZ, HDZ, MDZ, EIZ and UDZ.

Theoretical basis of the criterion for dividing EDZs based on R Ud
It is well known that the deviatoric stress-axial strain relationship curve of rock is based on a specific confining pressure. Thus, the resulting energy dissipation degree R Ud -axial strain relationship curve is also based on a specific surrounding pressure, while the surrounding rock's different position corresponds to different stress states. Therefore, the energy dissipation degree R Ud -axial strain relationship curve at a specific surrounding pressure should correspond to the energy dissipation state at any location in the surrounding rock. Fig 6 shows five curves of the relationship between the energy dissipation degree R Ud and the axial strain of granite in the Shuangjiangkou Hydropower Station at different confining pressures obtained by conventional triaxial tests. Through curve fitting, the five curves can be expressed with a simple monotonic increasing function, as shown in Eq (9).
The relevant fitting coefficients A 1 , A 2 , A 3 and A 4 and their values are shown in Table 2. Fig 6 shows that all the curves are well fitted, and the correlation coefficients R 2 reach above 0.96. As shown in Fig 7, the coefficients A 1 , A 2 , A 3 and A 4 are fitted to the confining pressure σ 3 in Eq (9) to represent the energy dissipation degree at different stress states in the surrounding rock. Eq (10) is the fitting equation.
where l i , m i and n i are the relevant fitting coefficients, and their values are shown in Table 3. Fig 7 shows that the correlation coefficients R 2 all reach above 0.98, which indicates that the coefficients A 1 , A 2 , A 3 and A 4 all have a clear parabolic-type relationship with the confining pressure σ 3 .
Eqs (9) and (10) are empirical formulas for the criterion for dividing the surrounding rock EDZs within the underground cavities based on the energy dissipation degree R Ud, derived from the numerical analysis of the results of conventional triaxial tests, which take into account both the effects of the surrounding pressure and strain, and they can be used to calculate the value of R Ud during rock damage at any position in the surrounding rock. Their formats are simple and can therefore be easily applied in engineering.
In practical engineering applications for a specific project, only a set of conventional triaxial tests needs to be carried out, and the empirical formulas for the criterion can be obtained by simple data processing according to the method provided in this paper. In particular, during the numerical simulation calculation, the stress-strain state of any unit in the model is known, and then the value of the energy dissipation degree R Ud of the unit can be calculated by Eq (9) and Eq (10). The dissipation energy in the numerical simulation can be calculated by Eq (8) and Up d = -0.044+0.015σ 3 (Fig 4(F)). The change in the surrounding rock's R Ud can be acquired, and the surrounding rock's damage degree and range can be divided through quantitative analysis of R Ud , which can provide reference advice on the design of support measures and then analyze and evaluate the surrounding rock's overall stability.

EDZs in the surrounding rock mass of underground caverns based on R Ud
The calculation object of this numerical simulation model is the left bank mountain where the underground workshop of Shuangjiangkou Hydropower Station is located, as shown in

PLOS ONE
The criterion for dividing the surrounding rock EDZs of underground caverns There are also three defects distributed in the area, including Fault 1, Fault 2 and a lamprophyre vein. As shown in Fig 8(B), the excavation of the underground workshop caverns is carried out in 9 stages, and the site construction is currently excavated to the 4th stage. The 3D model of the underground caverns is shown in Fig 8(C). The software used for the numerical simulation in this paper is FLAC 3D , and the yield criterion is the Morh-Coulomb criterion built into the software. The deformation modulus E d , bulk density γ, Poisson's ratio μ, cohesion c and internal friction coefficient f for numerical simulations are determined by the designer of the project, as listed in Table 4.
The inversion of the in situ stress field based on the measured in situ stress indicates that the in situ stress field in the engineering region is mainly composed of tectonic stress and gravitational pressure, and its distribution is influenced by the tectonic stress, geological structure,  and local topography [57]. From the top to the bottom of the surface, σ x , σ y and σ z increase with depth. After the fourth stage of excavation is completed, the minimum principal stress and maximum principal strain of all units in the model can be obtained and brought into Eqs (9) and (10) to derive the value of R Ud . Based on the four thresholds of R Ud in Table 1, the distribution ranges of the CDZ, HDZ, MDZ and EIZ of the underground caverns can be derived, as shown in Fig 9. The CDZ is mainly distributed in the area of 0.5~3 m around both the cavern and the connection area of the cavern. The HDZ is more widely distributed on the basis of the CDZ when the depth is increased to approximately 1~7 m, and the increase is more obvious in the sidewalls of cavities and the connection area between cavities. The MDZ is increased to approximately 1~12 m on the basis of the HDZ, especially the middle of the cavity sidewall, and the connection area between cavities is increased. The EIZ forms a large interconnected area at the periphery of the MDZ, and it is evident that the depth of the impact of the excavation on the surrounding rock is very large, even greater than 30 m at the deepest point.

EDZs of the surrounding rock mass of underground caverns based on wave velocity
The propagation velocity of elasticity waves varies in different media, and elasticity wave testing techniques have been widely utilized to study EDZs [24,26]. The rock integrity index (K v ) is usually defined as the square of the ratio of the elastic longitudinal velocity of the rock mass (v p ) to the elastic longitudinal velocity of the intact rock (v pr ), as shown in Eq (11). The higher the wave velocity of the longitudinal waves in the rock mass, the higher the value of K v , indicating the better quality of the rock. D is defined as the degree of damage to the rock mass, and the formula is shown in Eq (12); the higher the value of D is, the higher the degree of damage to the rock mass.
According to the variation curve of v p with borehole length, the EDZs of the surrounding rock can be divided into several zones, such as the CDZ, HDZ, MDZ, EIZ and UDZ, as shown

PLOS ONE
The criterion for dividing the surrounding rock EDZs of underground caverns in Fig 10. In the CDZ, v p is lower and becomes larger slowly with the borehole length, and D decreases slowly. The rock mass in the area is severely damaged and ruptured, and macroscopic fractures are developed. The rock mass loses bearing capacity or has very poor bearing capacity. The deformation of the surrounding rock is serious, and even rock explosion, flake gang, and crumbling occur. The v p is still low in the HDZ, but it increases rapidly with the borehole length, while D decreases rapidly. The rock damage and rupture in the area is more serious, the macroscopic fractures are more developed, and the bearing capacity of the rock is poor. In the MDZ, v p slowly increases to a peak with the borehole length, and D slowly decreases to a very low value. The rock mass in the area has a certain degree of damage rupture, and the majority of the fissures are not penetrated, while the rock has a certain bearing capacity. The v p and D variations in the EIZ and UDZ are not obvious and are difficult to use as a basis for the division of EDZs, so the EIZ and UDZ are unified into one region when the EDZs are divided based only on the wave velocity.
After the completion of the excavation of the fourth stage of the underground caverns of Shuangjiangkou Hydropower Station, eight wave velocity test boreholes were arranged in the 0+070.00 m section and 0+100.00 m section for the main powerhouse respectively. Based on the measured wave velocity test data, the EDZ distribution characteristics of the surrounding rock in the 0+070.00 m section and 0+100.00 m section of the underground caverns can be obtained as shown in Figs 11 and 12.
The distribution characteristics of the EDZs based on the elastic wave tests show that the depth variation range of the CDZ of the underground caverns at the 0+070.00 m section of the surrounding rock is 0.6~3.5 m, the depth variation range of the HDZ is 1.8~3.9 m, and the depth variation range of the MDZ is 3.7~5.1 m. The depth variation range of the CDZ of the underground caverns at the 0+100.00 m section of the surrounding rock is 1.0~2.3 m, the depth variation range of the HDZ is 1.5~3.3 m, and the depth variation range of the MDZ is 3.3~6.0 m. The wave velocity in the EIZ region and the UDZ region cannot be accurately determined with the variation in drilling depth, so the EIZ and UDZ can be combined into one region. Based on the analysis of the variation characteristics of wave velocity with drilling depth, the CDZ, HDZ and MDZ can be used as the key areas for support.  Table 5 shows the distribution characteristics of EDZs of the surrounding rock in 0 +070.00 m section of the main powerhouse based on R Ud thresholds and based on wave velocity. In the borehole ZKCD1-1 in the top arch area, the distances from the borehole entrance to the boundaries of the CDZ, HDZ and MDZ obtained with the R Ud -based method are 1.5 m,  Table 5 shows the distribution characteristics of EDZs of the surrounding rock in 0+100.00 m section of the main powerhouse based on R Ud thresholds and based on wave velocity. In the borehole ZKCD2-1 in the top arch area, the distances from the borehole entrance to the boundaries of the CDZ, HDZ and MDZ obtained with the R Ud -based method are 1.  The distribution characteristics of the EDZs of the surrounding rock of the underground caverns classified by R Ud -based and wave velocity-based tests show that the depths of the EDZs of the surrounding rock at the cavern sidewalls are greater than those of the top arch, the depths of the EDZs in the middle of the high sidewalls are significantly greater than that of the remaining parts, and the depths of the EDZs in the area affected by Fault 1 are greater than those in the unaffected area. The variation patterns of the EDZs of the surrounding rock classified based on wave velocity tests in the three major caverns in the 0+070.00 m section are generally similar to those based on R Ud . The distribution range of EDZs based on the wave velocity tests is generally slightly smaller than that of EDZs based on R Ud , which is related to the unfinished excavation of the caverns and the low number of wave velocity tests. The distribution characteristics of the underground cavern EDZs are obtained by numerical simulation calculations and in situ elastic wave tests, and the results are basically the same. However, the R Ud -based method for classifying the surrounding rock's EDZs has the obvious advantage of

PLOS ONE
The criterion for dividing the surrounding rock EDZs of underground caverns  being able to probe the boundaries of the surrounding rock's undamaged zone (UDZ) more explicitly, while the method based on wave velocity testing is not sufficiently explicit.

Discussion
The proposed energy calculation parameter R Ud represents the degree of dissipation energy change during the damage of the rock, yet there are some other published energy indexes for rock, such as the Local Energy Release Rate (LERR) which refers to the sudden release of energy stored in the rock mass per volume when the strain energy concentrated in a local area of the rock is larger than its limiting capacity during excavation. The formula for LERR can be written as [58,59]: where U emax and U emin are the peak and trough values of elastic strain energy intensity before and after brittle failure of the element respectively. LERR studies the release value of the elastic strain energy of the local surrounding rock after brittle damage, which is applied in numerical simulation calculations to predict the strength of the rock burst and the depth of the burst crater, while R Ud studies the damage zoning of the surrounding rock from the perspective of dissipation energy. The R Ud -based criterion can divide the EDZ into five zones, including CDZ, HDZ, MDZ, EIZ, and UDZ, while the range of rock bursts studied based on LERR is equivalent to the CDZ based on the R Ud -based criterion. In addition, the damage zoning of the surrounding rock based on R Ud can provide support design advice for the excavation of the

PLOS ONE
The criterion for dividing the surrounding rock EDZs of underground caverns surrounding rock, such as the support method, the length of the free section and anchor section of the prestressing anchor, etc.
Referring to the conventional triaxial test data for brittle sandstone [59], the threshold values of R Ud for different confining pressures were analyzed by applying the method provided in this paper and are shown in Fig 15 and Table 6. It is found that the same threshold value of R Ud for brittle sandstone under different confining pressures varies within an interval that can be approximated as a constant, but is different from that of granite. Therefore, in the application of the method of this paper, the threshold values of R Ud need to be determined on the basis of experiments for different rocks.
The criterion proposed in this paper is based on the analysis of conventional triaxial loading test data of intact granite, without considering the influence of unloading effect and existing fissures, which are important factors affecting the stability of underground cavern excavation [60][61][62], and we will make corresponding research in related fields in the next step.

Conclusion
The paper studied the energy dissipation mechanism and energy evolution characteristics during conventional triaxial tests of granite in Shuangjiangkou Hydropower Station, proposed the criterion for dividing the EDZs of the surrounding rock mass of underground caverns based on the energy dissipation degree R Ud , analyzed the stability of the surrounding rock in the underground cavern of Shuangjiangkou Hydropower Station, came to the following conclusions.
1. The four stress thresholds (σ cc , σ ci , σ cd and σ p ) correspond to the four R Ud thresholds (R Udc , R Udi , R Udd and R Udp ) on the energy dissipation degree-strain curve. The stress threshold at different confining pressures becomes larger with increasing confining pressure, while the corresponding threshold of energy dissipation degree R Ud varies in a smaller interval at different confining pressures. The same threshold has good consistency under different confining pressures, and R Udc can be taken as 0.02, R Udi as 0.17, R Udd as 0.34, while R Udp is 1. The criterion uses four energy dissipation degree thresholds (R Udc , R Udi , R Udd and R Udp ) to classify the surrounding rock mass into five categories: CDZ, HDZ, MDZ, EIZ and UDZ.
2. Based on the criterion, the EDZs of the surrounding rocks of the underground cavern group of the Shuangjiangkou Hydropower Station are analyzed. The CDZ is mainly distributed in the area of 0.5~3 m around both the cavern and the connection area of the cavern. The HDZ is more widely distributed on the basis of the CDZ when the depth is increased to approximately 1~7 m, and the increase is more obvious in the sidewalls of cavities and the connection area between cavities. The MDZ is increased to approximately 1~12 m on the basis of the HDZ, especially the middle of the cavity sidewall, and the connection area between cavities is increased. The EIZ forms a large interconnected area at the periphery of the MDZ, and it is evident that the depth of the impact of the excavation on the surrounding rock is very large, even greater than 30 m at the deepest point.
3. The distribution characteristics of the EDZs of the underground caverns of the Shuangjiangkou Hydropower Station at the 0+070.00 m section and the 0+100.00 m section of the surrounding rock are obtained based on in situ wave velocity tests and are compared with those based on R Ud . The R Ud -based method for classifying the EDZs of the surrounding rock has the obvious advantage of being able to probe the boundaries of the UDZ more explicitly, while the method based on wave velocity testing is not sufficiently explicit. Generally, the CDZ, HDZ and MDZ should be used as the focus areas for support.